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Abstract 

We estimate the effective heavy quark and antiquark potential in the quark gluon plasma using 
the gravity dual theory. Two models are considered: AdS§ and Sakai-Sugimoto model. The effective 
potential, obtained by using the rotating fundamental open-string configurations, has the angular 
momentum dependence which generalizes the static central potential. For zero angular momentum 
case, we obtain asymptotic form of the potential for general metric n, and for both gravity dual models 
where n = 3, 4. The mass dependence of the potential is derived at the leading order together with its 
temperature dependence. Motivated by the asymptotic form for zero angular momentum state, we fit 
the effective potential for J = 1, 2 states in the binding region. The fitting parameters are found to be 
functions of temperature. Finally, we discuss the differences and similarities of the effective potential 
between the two gravity dual models. An interesting result is that position of the minimum of the 
potential is determined only by angular momentum and independent of the temperature. 
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I. INTRODUCTION 



J/ip, a cc bound state, is a good probe for the formation as well as the properties of the 
quark gluon plasma (QGP). At low temperature, the interaction between c and c can be phe- 
nomenologically described as Coulomb-like potential 1/r. Because of confinement, one would 
expect there is additional linear potential ar which is dominant at long distance. When the 
temperature increases, the string tension a decreases. At the critical temperature where decon- 
finement phase transition takes place, one would expect that a ~ 0. In QGP, the existence of 
deconfined quarks and gluons would modify the interaction between quark and antiquark such 
that the static potential 

V{ T) = _^S2 e -r/r D P0 j (1) 

where r^, Debye screening radius, is the decreasing function of temperature. Roughly speaking, 
once the screening radius is smaller than the radius of J/ip, the quarkonium will dissociate. 
The dissociation of Jlip can be used as the signature of the formation of QGP as well as its 
temperature profile [1]. 

One can calculate the heavy quark and antiquark potential using the AdS/CFT correspon- 
dence . At zero temperature, the potential [4!] 

47T 2 v / A7 /l 



where A4 = gy M iV = 47rg s iV and r is the distance between quark and antiquark. One can see 
that the potential is proportional to —1/r. However, the confining potential ar is dominant 
at large distance in cold nuclear matter. Therefore we can see this model is more accurate 
to the deconfined phase where the a disappears, which requires the temperature to be higher 
than the deconfinement temperature T c . In order to take into account the temperature and 
thermal phenomena of the gauge plasma, one can consider a Schwarzschild-anti-de Sitter Type 
IIB supergravity compactification 

In the gravitational dual picture, the meson is obtained using the fundamental string whose 
two ends are connected with the probe brane. The two ends of this string represent the quark 
and antiquark separately. The position of the probe brane sitting at certain point in the radial 
direction is proportional to the mass of the quark. Since one expect that the presence of the 
probe brane would not affect the geometry in the leading order, it is necessary to put the 



probe brane far away from the branes which produce the gravitational background, and hence 
the large mass of the quark and antiquark. In the dual picture, there are two configurations 
which are competing with each other. One is the fundamental string connecting the quark 
and antiquark which corresponds to the bound state. The other is the two strings stretching 
from the probe brane to the horizon which represents the free quark and antiquark. At given 
temperature, one needs to compare the energy difference of these two configurations. It is used 
to determine the critical separation between quark pairs at which the bound states dissociate. 
The dissociation of bound state really implies the situation when the energy of the parallel 
strings is lower than the connected configuration. Then it is always preferred energetically 
for the pair to be separated as free quarks. Using this strategy, one would see that at finite 
temperature the potential between heavy quark and antiquark exhibits short-range asymptotic 
behaviour {f| which is qualitatively similar to the color-screened potential in Eqn. ([I]). 

The results from the gravity dual theory show that the static heavy quark potential has — 1/r 
behaviour at zero temperature and short-range asymptotic behaviour at the finite temperature. 
This results in scaling of the screening length L* with 1/T for states with zero angular momen- 
tum. It was also demonstrated that the screening length of the quark antiquark state scales 
with velocity of the moving meson as L* ~ (1 - v 2 ) l/A Q. However, one may ask if we can 
obtain the effective potential for the excited heavy quark bound states with nonzero angular 
momentum from the gravity dual theory. Excited heavy quark and antiquark bound states 
could have the non-zero angular momentum and spin. The potential between heavy quarks 



ar momentum, 
Therefore, 



which can move arbitrarily in the bound states would depend on the spin, angu 
as well as the velocity (e.g. see explicit form calculated using dual QCD in Ref. 
the potential would be more complicated than the static potential. Once the complete potential 
between heavy quark pairs is obtained, we can use it in the Schrodinger equation to get the 
complete spectrum of the mesonic states (see e.g. Ref. [3]). 

Due to the success of the gravity dual theory to describe the static potential of heavy 
quark pair, we are optimistic and hope to estimate the equivalent potential for the excited 
bound states from the gravity dual pic ture. To obtain the gravity dual picture, we boost 
the gravitational background metric [10] for two different models: AdS^ and Sakai-Sugimoto 
model. Taking advantage of the numerical technique, we calculate the effective potential of 
excited bound state as a function of the separation between quark pairs r, at fixed angular 
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momentum J and temperature T. The result shows that for each fixed angular momentum 
J and given temperature T, the effective potential becomes positive at small distance due to 
the angular momentum effect being dominant in this region. At the intermediate distance, 
the potential turns negative where the static binding potential dominates. Finally it becomes 
positive again at large distance because of screening effect. We also calculate the static potential 
for i = J = at nonzero temperature for two different models. The difference between the two 
models is discussed. 

This paper is organized as the following. In Section II, we demonstrate how introduction 
of the angular momentum barrier allows the effective potential to become positive at large 
distance, and therefore the existence of the "screening length". In Section III, we analytically 
and numerically obtain the effective potential for I = and excited £ = 1,2 bound states of 
heavy quark pair using the gravity dual models. The shape of effective potential produced 
by the gravity dual models are qualitatively similar to the expected 4-dimensional result from 
Section II. Significant differences between the two gravity dual models, n = 3,4, are discussed. 
In Section IV, fitting of the effective potential for excited £ = 1,2 states to the asymptotic form 
V e ff = C n — a n /r 2 /( n ~ 2 ) motivated from £ = case is performed and shows exceptionally good 
fit. We summarize our results in Section V. 

II. HEAVY QUARK AND ANTIQUARK POTENTIAL 

Because of the larger mass of heavy quark e.g. charm and bottom, the bound states of 
heavy quark and antiquark can be described using the potential model. One can expect that the 
heavy quark and antiquark bound state may be described by the solutions of the nonrelativistic 
Schrodinger equation 



where fi is the reduced mass of the quark- antiquark, with fi = m/2 for cc, bb quarkonium states. 
For the central potential, the wave function \l/(r) can be written as 
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FIG. 1: The effective potential between quark and antiquark in the quark gluon plasma for I = 
(left) and I = 1 (right) at finite temperature from Eqn. (JT]). Potential barrier induced by the 
orbital angular momentum can modify the effective potential so that the potential actually becomes 
positive at distance larger than a finite "screening length" L* as shown in the figure on the right. For 
demonstration purpose, we set rp = 10 for 1 = 1 figure. The screening length L* ~ 16 is substantially 
larger than the screening radius ro = 10 in this case. 

where R(r) is the radial wave function and Y( m is the spherical harmonic function. With this 
substitution, the radial wave function satisfies 

+ i)h 2 ' 



2 fx \dr 2 r dr J 



E-V(r) - 



2fir 2 



R(r) = 0. 



(5) 



One can define the effective potential 

V eff (r,£) = V(r) + 



2 fir 2 



(6) 



where I is the quantum number for the orbital angular momentum. The term induced by orbital 
angular momentum £(£ + l)h 2 /2fir 2 forms an angular momentum "barrier" ~ 1/r 2 dominating 
at small distance r. This angular momentum barrier is crucial when we consider fundamental 
differences between I = and I > quarkonium states when they are submerged into the QGP 
as we will see in the following. 

In the deconfined phase such as QGP, the linear confining potential or becomes zero. The 
remaining is the Coulomb-like potential which would be screened by the plasma. Analogous to 
the QED plasma, the screening QCD potential between quark and antiquark in the quarkonium 
submerged in the gauge plasma can be expressed as in Eqn. flTJ. The effective potential then 
becomes 

£{£+l)h 2 a(T) 



-r/r D (T) 



2fir 2 



(7) 



where rp is the Debye screening radius. 

The effective potential for £ = state approaches zero from below but never actually crosses 
over and becomes positive. On the other hand, the potential barrier induced when £ > 
modifies the total potential such that at finite distance L* > r^, the effective potential changes 
sign and becomes positive for r > L* as is shown in Fig. [TJ To distinguish L* where V e fj 
actually becomes zero from the Debye screening radius rp, we will call it "screening length" 
hereafter. 

The "melting" or dissociation of the bound state of quark and antiquark could be understood 
as the tunneling of the quark (antiquark) from the region where V e ff < 0, through the region 
of potential barrier V e ff > 0, to the large r region where it is effectively free. It is also possible 
that the thermal kinetic energy of the quark in the QGP (~ kgT) is substantially larger than 
the binding effective potential and thus resulting in almost complete dissociation of the £ > 
states (see e.g. Fig. [Q. This could explain why the £ = states are much harder to "melt" 
comparing to £ > states. 

We can see that the angular momentum induces significantly distinctive features to the 
effective potential and consequently to the behaviour of the excited quarkonium states. It is 
desirable to understand the more precise form of the potential as well as its dependence on the 
temperature so that we would gain more understanding of the properties of the QGP as well 
as the determining factors of the melting of quarkonia with varying angular momentum when 
they are submerged into the plasma. 

III. EFFECTIVE QUARK ANTIQUARK POTENTIAL FROM GRAVITY DUAL 
MODELS 

In Section II, one can see that nonrelativistic Schrodinger equation may be used to solve 
the heavy quark bound states provided that we have the correct form of the potential. The 
question is the exact nature of the screening potential between quark and antiquark in the 
deconfined phase at finite temperature. Analogy to QED plasma suggests the form of the 
screening potential in the form of the Debye screening, Eqn. ([!]). This form is expected from 



the perturbative calculation [12|] but there is no guarantee it would remain valid in the regime 
where the 't Hooft coupling is large. 

A complementary approach to calculate the effective potential when the gauge interaction 
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is strongly coupled with the large 't Hooft coupling is by means of the AdS/CFT correspon- 
dence |4| . One may ask if it is possible to obtain t dependent effective potential V e ff (r, £) in the 
gravitational dual picture. In the following part, we are going to obtain this kind of effective 
potential by boosting the background metric. 

In order to calculate the potential between the rotating quark and antiquark pair, we consider 
the following 5 dimensional metrics 

ds 2 = (^-Y /2 (-/„ (it) dt 2 + dp 2 + p 2 d<p 2 + dz 2 ) + ( ^) H/2 -pT-, (8) 

\R n J V ' \ U J Jn{u) 

where 

/r-M-l-i (9) 

u is the radial direction, z is the direction perpendicular to the plane of rotation and position 
of the horizon Uh = ^-^R^T 2 , tiR\T for n = 3, 4 respectively. 

The n = 3 metric, known as Sakai-Sugimoto model [3| (we ignore details of the compact 
5 dimensional manifold here assuming no Kaluza-Klein excitations with respect to these direc- 
tions) in the high temperature phase [lij], is the near-horizon limit of the metric induced by 
a configuration of N c D4-branes intersecting with Nf D8-branes and Nf anti-D8-branes. The 
(anti-)D8-branes are located at £4 = 0, (L) of the compactified 24 around which the D4-branes 
wrap. The string theory in this background is dual to the maximally supersymmetric SU(N C ) 
Yang-Mills theory in 1 + 4 dimensions with one dimension compactified on a circle with radius 
R 3 . 

The n = 4 metric is the near-horizon limit of the metric induced by a configuration of N 
D3-branes in the 1+9 dimensional background. In this limit, the non-compact 5 dimensional 
subspace is approximately AdS§. The string theory in this background is dual to the supersym- 
metric Af = 4 Yang-Mills in 1+3 dimensions when the other compact 5 dimensional manifold 
is S5. The horizon in the gravity picture induces the scale and Hawking temperature into the 
theory Q]. These quantities are identified as Aq C d and temperature of the quark gluon plasma 
in the dual gauge picture. 

In order to calculate the potential between quark and antiquark in the gravity dual picture, 
it is suffice to consider the classical action of the string configuration. 

The Nambu-Goto action is given by 

S = — ( dTdaJdet(G M Nd a X M daX N ). (10) 

Z7T J V 
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From the above background metric in Eqn. ©,©5 the generic action for a general worldsheet 
gauge, t = t, u = u(a), p = p(cr), f = ust, is 



S = — / dtda , 
2tt 



u 



12 



+ 



(£)>)t>»-^ 



11) 



\ \ fn (it 

where u', p' is the derivative with respect to a. 

In order to calculate the potential and solve the equation of motion, we assign the following 
ansatz for the rotating connected-string configuration, 



t = r, p = a, u — u(p), ip = u)t. 



The Nambu-Goto action becomes 



S = — dtdp, 

Z7T 



u' 2 f u 
t.(u) ViL 



(/n-^V), 



(12) 



(13) 



where u' is the derivative with respect to p. 

Then the potential for the connected string configuration whose ends located at p = po is 



V 



1 ff>o 



7T JO 



dp\ 



It 



12 



u 



u 

Rr, 



(14) 



The conserved angular momentum can be calculated from dL/du as 



j 



dp 



(15) 



We will identify the classical angular momentum J of the classical string with the orbital 



angular momentum quantum number i of the quarkonium (instead of + 1)) and use them 
interchangeably in this article. The quantum mechanical spin of the quark and antiquark in 
the quarkonium will be ignored since we consider only the bosonic string in the gravity dual 
picture. 

The corresponding regulating potential for the parallel strings configuration can be calculated 
in another worldsheet gauge, 



t = r, u = a, p = p 0} if — oot. 



The potential then becomes 



V ' = ~ X 1 7 

TT Ju c \ J, 



du 



u 



(16) 



(17) 



where u c = Uh/{1 — PqCU 2 ) 1 /™, the minimum distance in the u-direction for rotating parallel 
strings. It is interesting to note that in the rotating metric, the position of horizon where 
parallel strings end is shifted from Uh to u c = Uh/(l — plco 2 ) 1 ^ > Uh- The equation of motion 
in this gauge for parallel strings configuration is trivially satisfied since dL/d(d a u) = L for the 
Lagrangian L. 

On the other hand, the equation of motion for the connected-string configuration is 



d 
dp 



( \ 



f I ( u' 2 . ,2 n 2 , / u \ n \ , nu"- 1 f £ , ,2 n 2\ 



0. 



(18) 

For general situation where u ^ 0, we can numerically solve the equation of motion for a given 
boundary condition u(po) = u max ,u'(po) — > oo and given value of the ui. The relationship 
between V e ff and other quantities such as J and r = 2po can be obtained by eliminating the 
parameter u. 

In this article, we set u max = 20, R n = 1 throughout our numerical analysis. The u max (mass) 
dependence of V e ff for £ = case is given in subsection IIII Al The dependence of the binding 
part of the potential on R n , for u max — > oo limit, is R^, R\ for n = 3, 4 respectively. The choice 
of R n determines the strength of the binding potential of the quarkonium in the QGP and as 
we will see in Eqn. (1261) . Ra = 1 leads to ~ 0.23 (of the binding potential of the form — a^/r) 
which is close to the value used in conventional potential models [l|. 

A. The potential for lo = (£ = 0) case 



For this case, we have additional constant of the motion since the Lagrangian does not 
depend on p explicitly, 



,dL _ _ fn (£) 



du 

from which 



A = L-u'—= ; w (19) 

+ fn (it) 



U' 2 



n n 

Umax f?J 7/ 2 



r = 2 1 S^E du, (20) 
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where Ub = (u^ + Uq)" ,Uq = A 2 R™, and u max is position of the probe brane. The effective 
potential V — V, from Eqn. (fTi]) and Eqn. (IT7j) . can then be expressed as 

1 / /"Umax . /n n — 7/™ fllraoi \ 

V-V=i / ^=^A,-/ *, . (21) 

7T \y« 6 Y « — U b Ju h ) 

We will first consider the effective potential in u max —>■ oo limit and then calculate the leading 
order u max (mass) dependence subsequently. 
• Effective potential in u max — > oo limit: 



For T = in the u max — > oo limit of this u = (£ = 0) case, the regulated potential becomes 



1 r(l-i) \2JnT(l-h' 



(22) 



The potential has the form V ~ — A 5 /r 2 , — y/X~l/r with the 't Hooft coupling As 5 4 = 277r_R|/2, i?| 
for n = 3,4 respectively, in agreement with Ref. [l5l |. 

For nonzero temperature in M ma x - ► oo limit, by using asymptotic expansion as in Ref. 6j, 
we can obtain the leading order contribution with respect to A as (see generic form for general 
n in Appendix A) 

\l/(n-2) 

V eff (r, T) = C n (T) - a' n ^j^ [l + 0((rT) 2n ^)} , (23) 

where A = X^^ for n = 3,4 and C n (T) is a constant which is the increasing function with 
respect to the temperature and that C n (0) = 0, 

C n {T) = ^ (24) 

7T 

^ A 5 T 2 , ^aIt for n = 3, 4 (25) 



and 



:-!2 



v ^)(r(|)r(|)/r(|)) 3 ~ 0.00963 for n = 3, 

|(r(!)r(|)/r(i)) 2 ^ 0.2285 forn = 4. 

The other terms of the potential are of the positive power of rT and therefore suppressed as 
long as r < 1/T. From Eqn. (f23|) . it is obvious that V e ff(r = L*,T) = gives, at the leading 
order, 



L*T = ^(3/2)5, 0.2704, 0.2285 torn = 3, 4. (26) 
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FIG. 2: The effective potential for £ = J = as a function of the distance between quark and antiquark 
for the metric n = 3 (left) and n = 4 (right) for the £ = states. The solid line is for T = 0.20 GeV, 
the dashed line is for T = 0.25 GeV, and the dashed-dot line is for T = 0.30 GeV. 

Note that the screening length L* is independent of A regardless of n. 

The numerical plots of the effective potential at nonzero temperature for each gravity dual 
are presented in Fig. [2J We found that the fits to 1/r 2 , 1/r for n = 3,4 are valid with great 
accuracy for the physical range of r close to the screening length into the binding region. The 
asymptotic expansion works very well as long as the physical range of r under consideration is 
not too large. This provides a strong motivation to fit the form Eqn. (123]) to the binding region 
of more complicated cases with I > as we will see in Section IV. 

• u max (mass) dependence of the effective potential: 



We have set u max — > oo for the above analytical calculations. The limit u max — > oo provides 
universal form of potential for heavy quarkonia in the very large mass limit. On the other hand, 
the specific mass dependence of the potential can be achieved when we fix u max (~ m, mass of 
the quark) to some large finite value. 

From Eqn. (J2T1) . we can explicitly express M maa; -dependent piece of V e ff as 

Veff = V eff (u max — > oo) + V eff (u max ) (27) 
1 / r°° y/u n - ul 

7T \JUmax \/U n Ufo J Umax 



V ef f(u max ) = --[ I lnn _ % du- I du) (28 




\A - (WV) 



n POO \ 

= dy- J dy\ , 



where h,b = Uh,b/u r , 



(29) 
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By expanding with A = h,b } 
(1 - (A/yD^ 2 

we obtain 

We can eliminate b by using Eqn. (120]) , 



Rn Uq 



it 
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- L 2V- 1 -2 
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■ yu lb j 
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,h n b n 



2n(n - 1) An(2n - 1) 



+ ... . 
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du 



U — M 



[U 



' 1— n 



2 J R n / 2 M^- n )/V6" - 



-l/n 



O E.n/2 (2-n)/: 



b (2-n)/2 



n 



A 



1 



6 



O n 2-^1 



n 



1 1 3 
~'2'2 



1 

n 



+ 0(m 



1— n \ 
max) 



(30) 
(31) 

(32) 
(33) 

.(34) 



where a n is given in Appendix A. Taking the leading contribution with respect to 0((h/by 
and inverting to obtain 



1 



2ar,R 



n/2\ 2n/(n-2) 



-/i n [1 + 0{(h/b) n )} + .. 



.. . , ,. ... . (35) 

(n-2)(3n-2) L vw y /J v ; 

We substitute Eqn. fl35l) into Eqn. fl3~TT) to obtain the leading order (with respect to u max ,u h ) 

dependence of effective potential on the mass (defined at T = 0), m = u max /2-n as 



Veff{m) 



m 



l-n 



in 



4(n 



-u h - 



2a n R 



t/2 \ 2n/(n-2)- 



nr 



(36) 



l)(27r)» Un-2)(3n-2) h 

The mass dependence contribution is sufficiently suppressed even at u max = 20. For n = 3, 4, the 
leading order mass dependence V e ff(m) ~ 0(l/m 2 ), 0(l/m 3 ) respectively. The r-dependence of 
V e ff(m) is interestingly of order 0(l/r 6 ), 0(l/r 4 ) and the temperature dependence is ~ x 6 ,^ 
for n = 3, 4 metric model. It is interesting to compare this result with the mass-dependent 
potential calculated in other approaches such as pNRQCD (see Ref. 16) and references therein). 
Note that our result, Eqn. fl36l) . is applicable when the quarkonium is submerged in the QGP. 

For (u) I > case, since we are more familiar with the Coulomb-like potential ~ 1/r from 
the n = 4 AdS§ gravity dual, we will first consider the effective potential calculated in this case 
in subsection IIII Bl The effective potential calculated in the n = 3 Sakai-Sugimoto model will 
then be considered subsequently in subsection IIII CI Comparison between the two gravity dual 
models will be discussed in subsection IIII Dl 
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B. Effective Potential from n = 4 metric, AdS§ Model 



From Eqn. (|22|) . we can derive the well known result for the Coulomb-like potential at 
zero temperature for £ = state, V e ff ~ — V%t/ r - At nonzero temperature, the asymptotic 
expansion also gives V e // ~ — V%i/ r as the leading order contribution with respect to the large 
't Hooft coupling A4 [6|. The next-to-leading order C&(T) term is a constant proportional to the 
temperature T {uh(n = 4) ~ T) and the other terms are of positive power of rT and therefore 
suppressed as long as r < 1/T. 

As in Eqn. ((7j), we expect that inclusion of the angular momentum barrier brings in pos- 
itive 1/r 2 term dominant at small r as well as reducing the value of the screening length L* 



when compared with £ = case at the same temperature (see Table I of Ref. [17|). The ef- 
fective potential of the states with I = J — 1, 2 at varying temperature is given in Fig. [9], [TUl 
respectively. 

At T — 0, £ — 1, the potential starts with positive value at small r due to the angular 
momentum effect and turns negative around r = 0.04 fm. The effective potential becomes 
minimal at r = 0.06 fm and becomes less negative at larger distance but never becomes positive. 
It is obvious that as temperature gets higher, the potential becomes weaker and the screening 
length gets shorter. Interestingly, the distance tq where the potential is minimum is the same 
regardless of the temperature. Position of the minimum, r , is determined only by the angular 
momentum J of the string. As we will see later on, this interesting behaviour also exists in the 
case of n = 3 Sakai-Sugimoto model for both £=1,2 cases (as long as T > T c ). 

The behaviour of £ = 2 is very similar but somewhat fading out comparing to £ = 1 case. 



This is the sign that the state with £ = 2 is less tight 
than £ = 1. We can see that from Fig. 5 of Ref. 



y bound and it is much easier to dissociate 



£ = 2 state melts around T = 0.24 GeV 



while £ = 1 does not melt as T goes as high as 0.45 GeV or so (we will see later in Fig. [3] that 
it melts around T = 0.59 GeV) 1 . It should be emphasized that the position of minimum r is 
fixed and independent of temperature as in £ = 1 case. 



1 Even though the effective potential is negative in certain region, it is not necessarily true that it always admits 
quantum mechanical bound states in that region (e.g. if the binding region is too shallow). Therefore quarko- 
nium, in practice, start to dissociate at lower temperature than the melting temperature determined classically 
when V e ff > (for all r) as in Fig. [3) We can think of the classically determined melting temperatures as 
the upper bounds on the melting temperatures above which complete dissociation is guaranteed. 
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The numerical results show that for given temperature T, one encounters two different 
regions for r > tq relevant to the melting of quarkonium. When the separation between quark 
and antiquark r < aj/T(where aj is the constant determining the screening length depending 
on J), there is the Coulomb-like behavior. When r > aj/T, one can see that the potential 
would cross over the r-axis and becomes positive. The quark and antiquark become free and 
the bound state dissociates. 

C. Effective Potential from n = 3 metric, Sakai-Sugimoto Model 

At zero temperature for £ = state, the potential from this metric is V e ff ~ — As/r 2 . Even 
though this is not necessarily physical since we expect a phase transition to different metric 
when the temperature T drops below T c , this form of potential provides hint to what kind of r 
dependence the effective potential has for generic £ at nonzero temperature. 

At nonzero temperature, the asymptotic expansion for £ = gives V e ff ~ — A/r 2 + C$(T). 
The next-to-leading-order term C3 is proportional to T 2 (iih(n = 3) ~ T 2 ) as opposed to T of 
the n = 4 case. 

For £ > 0, we present the results in Fig. [TTJ [12] for £ = J — 1, 2 respectively. The angular 
momentum barrier is dominant at small r. The effective potential turns negative around r = 
0.05 ~ 0.06 fm for £ = 1, and around r = 0.085 ~ 0.095 fm for £ = 2. The minimum of potential 
is at r = 0.065 (0.11) fm for I = 1 (2). The melting of £ = 1, 2 occurs at T = 0.50, 0.35 GeV 
respectively as we can see in Fig. [31 

The general shape of the potential curve for £ = 2 does not change much from £ = 1 in 
this metric model in contrast to the case of n = 4 metric. This is due to the smaller curvature 
with respect to the radial direction u of Sakai-Sugimoto model. The state with higher angular 
momentum will be rotating at larger distance in u coordinate and the effect of curvature will 
be seen more distinctively between the two gravity dual models. We can see that the melting 
temperature for 1 = 1 states in n = 3 model is around 0.35 GeV, considerably higher than that 
of n = 4, around 0.24 GeV, as is shown in Fig. [3j The I = 2 states of n = 4 metric model 
dissociate much earlier and easier. 

From Fig. [TTlfT2| we can see that the screening length of this metric model is interestingly 
shorter than that of n = 4 model for £ = 1, and larger than the value of n = 4 for £ = 2 at the 
same temperature. 
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D. Comparison between the Two Metric Models 



The curvature of the dual metric relevant to the physics of the gauge plasma is the one with 
respect to the radial direction u. In our setup, the higher n naturally implies higher curvature 
and we expect the n = 4, AdS$ model to show certain enhancing effects when compared to 
the n = 3 Sakai-Sugimoto model. The rotating string with higher angular momentum will be 
spinning at further distance from the horizon and we expect to see effects of curvature more 
distinctive than states with lower angular momentum. 
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FIG. 3: The effective potential when the quarkonium melts, n = 3, £ = J = 1,T = 0.50 GeV and 
£ = J = 2, T = 0.35 GeV on the left and n = 4, £ = J = 1, T = 0.59 GeV and £ = J = 2, T = 0.24 
GeV on the right. 



These curvature effects can be seen in Fig. |3]|5l Fig. [3] shows the melting temperatures for 
£ = 1,2 states for each gravity dual model. While the melting temperature for I = 1 state of 
n = 4 is higher than that of n = 3, the melting temperature for £ = 2 state becomes drastically 
lower than that of n = 3 model. 
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FIG. 4: T = 0.20 GeV, £ = J = 1 (2) for n = 3, 4 on the left (right). 
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n 


J = 1 


2 


3 


0.065 fm 


0.11 fm 


4 


0.06 fm 


0.15 fm 



TABLE I: The values of distance ro where the effective potential V e ff(r,T) is minimum for the dual 
metric n = 3,4. 



In Fig. HI it is obvious that n = 3 potential binds stronger than the potential of n = 4 metric 
model, with much higher binding energies. We can see that even at T = 0.20 GeV, the £ = 2 
state of n = 4 already almost melts away while the same state of n = 3 is still strongly bound. 




0.1 0.15 0.2 0.25 u.uj u.-l u.-lj 

r(fm) r(fm) 

FIG. 5: T = 0.20 GeV, £ = J = 1, 2 for n = 3, (4) on the left (right). 



Fig. [5] represents comparison between states with i = 1,2 for each metric. The separation 
of ro in n = 3 model is about 0.045 fm, much lower than in n = 4 where the gap of ro between 
I = 1, 2 is around 0.09 fm. The values of the minimum of potential r are summarized in 
TableE 

It is observed in subsection IIII CI that the screening length for I = 1 (2) from n = 3 model is 
smaller (larger) than the corresponding value from n = 4 model. Interestingly, from Eqn. fl26l 
and Fig. [21 the screening length of i = state from n = 3 metric model is also found to be 
larger than the value from n = 4 model. 

Finally, we also found numerically that the angular momentum J does NOT depend on the 
temperature T, only depends on u,r. We can see from Fig. l9lfT2l that for fixed value of J and 
u, r is determined to be the same for each temperature curve (points with the same uj are on 
the same r once J is fixed). It appears that submerging quarkonium into the QGP at varying 
T does not affect its angular momentum. In other words, the interaction between QGP and 
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quark (antiquark) is radial and preserving angular momentum. 

IV. FITTING OF THE EFFECTIVE POTENTIALS IN THE BINDING REGION 

From Eqn. (1231) in the case of t = 0, it is natural to wonder if this form of effective potential 
also describes the binding region of the more general i > cases. By considering the shape 
of V e ff in both metric models, a simple guess is that far away from the angular momentum 
barrier ~ 1/r 2 , the binding region < r < L* should be approximated to a good precision by 
the asymptotic form V e ff = C n (T) — a n (T) / V 2 /( n-2 ). We actually found that this is the case. 

».5[ 1 i ] o.5f — r — ~" — ~ — "~ ^ ^ H 




r(fm) r(fm) 

FIG. 6: The fitting of the physical binding region r$ < r < L* with V e ff(r,T) = C3 — a-^/r 2 for 
t = J = 1, n = 3 dual metric. 
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FIG. 7: The fitting of the physical binding region tq < r < L* with V e ff(r,T) = C4 — a^/r for 
t = J = 1, n = A dual metric. 



In Fig.[6]and[7l the fits to £ = J = 1 for both metric models are demonstrated at T = 0, 0.20 
GeV. We also found that the same fitting works very well for the case i = 2. Remarkably, the 
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fitting works well at any temperature for both £=1,2 for both metric models. The values of 
the best fit parameters for £ = 1, 2 are summarized in Table HTl II III respectively. 



1 j Lrev 


as/ (GeVJ(tm) 


CK4/ (vjeVJ(tmj 


C3/ vjev 




o 


0.00695 


0.04004 








0.10 


0.00709 


0.04126 


r\ r\ A on 

0.0480 


d 1 1 on 

0.11230 


0.20 


0.00648 


0.04037 


0.1591 


0.19646 


0.25 


0.00618 


0.0390 


0.2428 


0.23014 


0.30 


0.00596 


0.03758 


0.3484 


0.26211 


0.35 


0.00538 


0.03503 


0.4315 


0.28031 


0.40 


0.00447 


0.03194 


0.4802 


0.29134 



TABLE II: The values of the fitting parameters a n , C n for the binding region ro < r < L* for the 
dual metric n = 3, 4 for I = J = 1. 



T/ GeV 


o 3 / (GeV)(fm) 2 


o 4 / (GeV)(fm) 


C 3 / GeV 


C 4 / GeV 





0.01149 


0.04266 








0.10 


0.01135 


0.03779 


0.0450 


0.07936 


0.17 


0.01004 


0.03132 


0.1149 


0.11608 


0.20 


0.00926 


0.02639 


0.1465 


0.11995 


0.25 


0.00830 


N/A 


0.2226 


N/A 


0.30 


0.00623 


N/A 


0.2631 


N/A 



TABLE III: The values of the fitting parameters a n , C n for the binding region tq < r < L* for the 
dual metric n = 3, 4 for £ = J = 2. 

A few observations can be made regarding the fit values of a n (T),C n (T). For both I = 
1,2 cases, C n is an increasing function of the temperature T for both metric models. The 
dependence of a n on the T is a bit more complicated. For I = 2 in both models, a n (T) is a 
decreasing function of the temperature. On the other hand for i = 1, a n (T) increases with 
T until around T = 0.10 GeV, then drops steadily as T gets larger. Since the deconfinement 
temperature T c is larger than 0.10 GeV, we can say that in QGP phase, a n (T), for each £, n, is 
a decreasing function of T. This is in contrast to £ = case where a n is constant (containing 
't Hooft coupling). 
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FIG. 8: The temperature dependence of the fitting parameter C n (T) for I = J = 1 (left), 2 (right). 

In Fig. [HI the nature of increasing function of temperature C n (T) for i — 1,2 appears 
consistent with the fact that for i = 0, C^(T) ~ T 2 , T respectively. There is small deviation 
at large T from the I = form due to rotation effects. 

V. CONCLUSIONS AND DISCUSSIONS 

We have considered the effective potential between heavy quark and antiquark being sub- 
merged in the QGP for the states with I = 0, 1,2. We have shown that even in conventional 
color-screened potential, introduction of the angular momentum barrier makes the excited states 
of quarkonium much less tightly bound, as well as making the potential at nonzero temperature 
crossing over zero value at finite distances. We define the distance where the effective potential 
turns from negative to positive as distance grows as "screening length" L*, in contrast to the 
screening radius where the potential is suppressed exponentially but not exactly zero. 

The effective potentials for I = at zero and nonzero temperature at the leading order are 
derived analytically for general n (= 3,4) in our setup, Eqn. (1221) . (1231) (see also Appendix A). 
The asymptotic expansion works well as long as r < 1/T. The mass dependence of the effective 
potential, Eqn. (136]) . is derived at the leading order together with its temperature dependence. 
The leading order mass dependence is of 0(l/m 2 ),0(l/m 3 ) for n = 3,4 respectively. This 
dependence remains even at zero temperature and depends on r as 0(l/r 6 ), 0(l/r 4 ) for n = 3, 4. 

The numerical plots of the effective potential V e ff(r, T) are given for each metric model with 
£ = J — 0, 1, 2 at various temperature values. An intriguing observation is that the position of 
the minimum r of the effective potential is fixed once we fix the angular momentum, indepen- 
dent of the temperature. Only angular momentum determines r of the effective potential for 
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each metric model. 

Due to higher curvature, the MS* metric model giws "early-melting potentia. for states 
with high angular momentum. As is shown in Ref. [17[, while £ = 2 state of n = 3 model 
resists to melting until T = 0.35 GeV, £ = 2 of n = 4 melts as early as T = 0.24 GeV. The 
shape of V e ff for £ = 2 of n = 4 metric model also shows very lightly bound potential at low 
temperature. 

The angular momentum of the quark antiquark system shows independence with respect 
to the temperature of the QGP, depending only on u and r. This can be understood that as 
long as the interaction between quark and antiquark calculated from the gravity dual models is 
radial, angular momentum of the bound state will not be altered as the temperature changes. 

Finally, fitting of the asymptotic form V e ff = C n (T) — a n (T) /r 2 ^ n ~ 2 ^ (motivated from £ = 
case) to the potential of £ = 1, 2 states in both metric models works very well for zero and 
nonzero temperature. While the constant C n is an increasing function of temperature, the 
constant a n is a decreasing function for T > T c . 
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APPENDIX A: ASYMPTOTIC FORM OF EFFECTIVE POTENTIAL FOR £ = 

Since n could be related to p of the Dp-branes which act as the source generating the curved 
background metric of the gravity dual as n = 7 — p, it is useful to express the asymptotic form 
of the effective potential V e ff(r,T) for £ = in u max — > oo limit, Eqn. (123]) . as function of 



generic n, 




(Al) 



where 




(A2) 
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FIG. 9: n = A,£ = J = l,T = 0, 0.20, 0.25, 0.30, 0.40 GeV from bottom to top respectively. 
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FIG. 10: n = 4, i = J = 2, T = 0, 0.17, 0.20, 0.30 GeV from bottom to top respectively. 
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FIG. 12: n = 3,1 = J = 2,T = 0, 0.17, 0.20, 0.25, 0.30 GeV from bottom to top respectively. 
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